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Abstract 

We consider a pure 2m-qubit initial state to evolve under a particular quantum me- 
chanical spin Hamiltonian, which can be written in terms of the adjacency matrix of 
the Johnson network J(2m,m). Then, by using some techniques such as spectral dis- 
tribution and stratification associated with the graphs, employed in [1, 2], a maximally 
entangled GHZ state is generated between the antipodes of the network. In fact, an 
explicit formula is given for the suitable coupling strengths of the hamiltonian, so that 
a maximally entangled state can be generated between antipodes of the network. By 
using some known multipartite entanglement measures, the amount of the entanglement 
of the final evolved state is calculated, and finally two examples of four qubit and six 
qubit states are considered in details. 

Keywords: maximal entanglement , GHZ states, Johnson network, Strat- 
ification, Spectral distribution 

PACs Index: 01.55.+b, 02.10.Yn 
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The idea to use quantum spin chains for short distance quantum communication was put 
forward by Bose [3]. After the work of Bose, the use of spin chains [4]-[18] and harmonic 
chains [19] as quantum wires have been proposed. In the previous work [1], the so called 
distance-regular graphs have been considered as spin networks (in the sense that with each 
vertex of a distance-regular graph, a qubit or a spin was associated) and perfect state transfer 
(PST) of a single qubit state over antipodes of these networks has been investigated. In that 
work, a procedure for finding suitable coupling constants in some particular spin Hamiltonians 
has been given so that perfect and optimal transfer of a quantum state between antipodes of 
the corresponding networks can be achieved, respectively. 

Entanglement is one of the other important tasks in quantum communication. Quantum 
entanglement in spin systems is an extensively-studied field in recent years [20,21,22,23], in the 
advent of growing realization that entanglement can be a resource for quantum information 
processing. Within this general field, entanglement of spinl/2 degrees of freedom, qubits, has 
been in focus for an obvious reason of their paramount importance for quantum computers, 
not to mention their well-known applicability in various condensed-matter systems, optics and 
other branches of physics. In [24], authors attempted to generate a Bell state between distant 
vertices in a permanently coupled spin network interacting via invariant stratification graphs 
(ISGs). At the first step, they established an upper bound over achievable entanglement 
between the reference site and the other vertices. Due to this upper bound they found that 
creation of a Bell state between the reference site and a vertex is possible if the stratum of 
that vertex is a single element, e.g. antipodal ISGs. The present work focuses on the provision 
of GHZ state, by using a 2m-qubit initial product state. To this end, we will consider the 
Johnson networks J(2m, m) (which are distance-regular) as spin networks. Then, we use the 
algebraic properties of these networks in order to find suitable coupling constants in some 
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particular spin Hamiltonians so that 2m-qubit GHZ state can be achieved. 

The organization of the paper is as follows: In section 2, we review some preliminary 
facts about graphs and their adjacency matrices, spectral distribution associated with them; 
In particular, some properties of the networks derived from symmetric group Sn called also 
Johnson networks are reviewed. Section 3 is devoted to 2m-qubit GHZ state provision by 
using algebraic properties of Johnson network J(2m,m), where a method for finding suitable 
coupling constants in particular spin Hamiltonians is given so that maximal entanglement in 
final state is possible. The paper is ended with a brief conclusion and two appendices. 

2 preliminaries 

2.1 Graphs and their adjacency matrices 

A graph is a pair F = {V,E), where y is a non-empty set called the vertex set and E is a, 
subset of {{x,y) : x,y & V,x ^ y} called the edge set of the graph. Two vertices x,y & V are 
called adjacent if {x,y) e E, and in that case we write x ~ y. For a graph F = {V,E), the 
adjacency matrix A is defined as 



otherwise 

Conversely, for a non-empty set V, a graph structure is uniquely determined by such a matrix 
indexed by V. The degree or valency of a vertex x eV is defined by 



where, | • | denotes the cardinality. The graph is called regular if the degree of all of the vertices 
be the same. In this paper, we will assume that graphs under discussion are regular. A finite 
sequence Xo,Xi, ...,Xn & V is called a walk of length n (or of n steps) if ~ Xi for all 
i — 1,2, Let l'^{V) denote the Hilbert space of C-valued square-summable functions on 




(2-1) 



i^i^) ^ \{y &V : y x} 



(2-2) 
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V. With each 13 & V we associate a vector \(3) such that the /3-th entry of it is 1 and all of 
the other entries of it are zero. Then {\/3) : /3 e V} becomes a complete orthonormal basis of 
1^{V). The adjacency matrix is considered as an operator acting in l'^{V) in such a way that 



2.2 Spectral distribution associated with the graphs 

Now, we recall some preliminary facts about spectral techniques used in the paper, where more 
details have been given in Refs. [26,27,28,29] 

Actually the spectral analysis of operators is an important issue in quantum mechan- 
ics, operator theory and mathematical physics [30,31]. As an example ii[dx) = \%l){x)\'^dx 
{lJ>{dp) — \ip{p)\'^dp) is a spectral distribution which is assigned to the position (momentum) 
operator X{P). The mathematical techniques such as Hilbert space of the stratification and 
spectral techniques have been employed in [32,33] for investigating continuous time quantum 
walk on graphs. Moreover, in general quasi-distributions are the assigned spectral distribu- 
tions of two hermitian non-commuting operators with a prescribed ordering. For example the 
Wigner distribution in phase space is the assigned spectral distribution for two non-commuting 
operators X (shift operator) and P (momentum operator) with Wyle-ordering among them 
[34, 35]. It is well known that, for any pair {A, \4>o)) of a matrix A and a vector \4>q), one can 
assign a measure /i as follows 



where E{x) — J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e.. 




(2-3) 



a~/3 



fi{x) = (0o|^(a:)|0o), 



(2-4) 




(2-5) 



Then, for any polynomial P{A) we have 




(2-6) 
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where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-4) and (2-6), the expectation value of powers of adjacency matrix A over 
reference vector |0o) can be written as 

{(l)o\A^\(l)o) = [ x'^liidx), 0,1,2,.... (2-7) 

JR 

Obviously, the relation (2-7) implies an isomorphism from the Hilbert space of the stratification 
onto the closed linear span of the orthogonal polynomials with respect to the measure /x. 



2.3 Underlying networks derived from symmetric group S,, 



n 



Let A = (Ai, A^) be a partition of n, i.e., Ai -I- ... -\- \m — n. We consider the subgroup 

Sm ® Sji-m of -S'^ with m < \^]. Then we assume the finite set = with IM^I = 

m!(n-m)! vcrtcx Set. In fact, is the set of (m — l)-faces of (n — l)-simplex (recall that, 
the graph of an {n — l)-simplex is the complete graph with n vertices denoted by Kn). If we 
denote the vertex i by m-tuple {ii,i2, ■■■,im), then the adjacency matrices A^, k — 0,1, ...,m 
are defined as 



2J 



1 if d{i,i)=k, 

, k^O,l,...,m. (2-8) 

otherwise (i,j e M-^) 



where, we mean by d{i,j) the number of components that i — (ii, ...,im) and j — (j'l, ...,jm) 
are different (this is the same as Hamming distance which is defined in coding theory). The 
network with adjacency matrices defined by (2-8) is known also as the Johnson network J{n, m) 
and has m -|- 1 strata such that 



Ko = 1, Ki 



m 



y m — I j 



I \ 

n — m 



, ; = l,2,...,m. (2-9) 



One should notice that for the purpose of maximal entanglement provision, we must have 
Km — 1 which is fulfilled if n = 2m, so we will consider the network J(2m, m) (hereafter we 
will take n — 2m so that we have Km = !)■ If we stratify the network J{2m, m) with respect to 
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a given reference node \4>q) = i2, ...,im), where i2, im) = |0...0^_1__^0...0^_1_^0...0^_1_^0) 
and ■■■ 7^ The unit vectors i — 1, ...,m are defined as 

\(p2) ^ —j= Y Yl \ii,...ii-i,i'i,ii+i,...,ik-i,i'k,ik+i---,im), 



' Kyi 



■3 kl^k2j^...^kj = l 

10^) = -^^ E (2-10) 

Since the network J(2m, m) is distance-regular, the above stratification is independent of the 
choice of reference node. The intersection array of the network is given by 

bi^{m-lf ; ci^P. (2-11) 

Then, by using the Eq. (B-49), the QD parameters ccj and Ui are obtained as follows 

ai^2l{m-l), Z = 0, l,...,m; uji^f{m-l + iy, Z = l,2,...,m. (2-12) 

Then, one can show that [27] 

A\(f>i) ^{l + l){m- + 2l{m - l)\(t>i) + l{m - I + (2-13) 

3 GHZ state generation by using quantum mechanical 
Hamiltonian in the network J(2m, m) 

The model we consider is the distance-regular Johnson network J(2m, m) consisting oi N — 
= lll^ sites labeled by {1, 2, N} and diameter m. Then, we stratify the network with 
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respect to a chosen reference site, say 1 (the discussion about stratification has been given in 
appendix A; In these particular networks, the first and the last strata possess only one node, 
i.e., |0o) = |1) and \4>rn) = 1^))- At time i = 0, a 2m-qubit state is prepared in the first 
(reference) site of the network. We wish to provide a maximal quantum entanglement between 
the state of this site and the state of the A'"-th site after a well-defined period of time, in which 
the corresponding network is evolved under a particular Hamiltonian. 

If the network be assumed as a spin network, in which a spin- 1/2 particle is attached to 
each vertex (node) of the network, the Hilbert space associated with the network is given by 
H = (C'2^®2m_ rjj^g standard basis for an individual qubit is chosen to be |0) = | J,), |1) = | t)- 
Then we consider the Hamiltonian 

Hs^l E (3-14) 

l<i<j<2m 

where, Hij — ui ■ Uj and Uj is a vector with familiar Pauli matrices erf, erf and erf. One can 
easily see that, the Hamiltoniaan (3.14) commutes with the total z component of the spin, i.e., 
[c^totaii-^s] — 0, hence the Hilbert space % decompose into invariant subspaces, each of which 
is a distinct eigenspace of the operator crf^^^;. So the total number of up and down spins are 
invariant under action of Hamiltonian or time evolution operator. Now, we recall that the kets 
12-1 ■■■ I i2m) with ii, . . . , i2m G {t) i} form an orthonormal basis for Hilbert space H. Then, 
one can easily obtain 

i/,,|...^...J^...) = |...J^...^...) 

i j i j 

and 

i/,,|...^...^...) = -|...J_^...^...) + 2|...^...^...). (3-15) 

i j i j i j 

Equation (3.15) implies that the action of on the basis vectors is equivalent to the action 
of the operator 2Pij — 7, i.e. we have 

Hij = 2Pij - I (3-16) 
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where Pij is the permutation operator acting on sites i and j. So 

/, (3-17) 



l<i<j<2m l<i<j<2m ^ 



^ 2m ^ 



V 2 / 



In fact restriction of the operator Yli\<i<i<2m Pij ^^e m-particle subspace (subspace spanned 

by the states with m spin up) which has dimension C^, is written as the adjacency matrix A 
of the Johnson network J{2m,'m), as 

Pij = A + m{m - 1)1. (3-18) 

l<i<j<2m 

For more details sec Rcf. [1]. Then we stratify the network with respect to a chosen refer- 
ence site, say \4>o). At time t — 0, the state is prepared in the 2m-qubit state \ip{t — 0)) = 
1 1 1 ... 1) I 00 ... 0) . Now, we consider the dynamics of the system to be governed by the Hamil- 

m m 

tonian 

m 

H^YJkPk{l/2 Y ^i-(^j + l^I)^ (3-19) 

k = l<i<j<2m 

Then, by using (3.17)-(3.20), the Hamiltonian can be written as 

m 

H^YWA) (3-20) 

it=0 

Jfe is the couphng strength between the reference site |0o) and all of the sites belonging to the 
k-th stratum with respect to \4>o), and Pk{A) are polynomials in terms of adjacency matrix 
of the Johnson network. Then, the total system is evolved under unitary evolution operator 
U{t) — e~'^* for a fixed time interval, say t. The final state becomes 

N 

m))-j:Mm (3-21) 
j=i 

where, N is the number of vertices, \j)s have 2m entries inclusive m entries equal to 1 and the 
other entries are and \A) = | 1^^^00_^^) so that /^^(t) := {j\e-'"^\A). 

m m 

The evolution with the adjacency matrix H = A = Ai lor distance- regular networks 
(see Appendix B) starting in \(f)o), always remains in the stratification space. For distance- 
regular network J(2m,m) for which the last stratum, i.e., contains only one site, then 
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maximal entanglement between the starting site \4>q) = 1^4) and the last stratum (the 
corresponding antipodal node) is generated, by choosing suitable coupling constants J^. In 
fact, for the purpose of a maximally entangled GHZ state generation between the first and 
the last stratum of the network, we impose the constraints that the amplitudes (0i|e~*^*|0o) 
be zero for all i — l,...,m — 1 and (0o|e~'^*|0o) — /, (0m|e~*''^*|0o) — f- Therefore, these 
amplitudes must be evaluated. To do so, we use the stratification and spectral distribution 
associated with the network J(2m, m) to write 



Let the spectral distribution of the graph is ii{x) — J2k=olkS{^ ~ ^k) (see Eq. (B-53)). 
The Johnson network is a kind of network with a highly regular structure that has a nice 
algebraic description; For example, the eigenvalues of this network can be computed exactly 
(see for example the notes by Chris Godsil on association schemes [39] for the details of this 
calculation). Indeed, the eigenvalues of the adjacency matrix of the network J{2m,m) (that 
is Xfc's in iJ,{x)) are given by 

Xk^m'^ -k{2m+l-k), A; = 0, 1, . . . , m. (3-22) 

Now, from the fact that for distance- regular graphs we have Ai — y/KlPi{A) [27], (0j|e~*^*|0o) — 
implies that 

m 

Y,lkPr{.Xk)e-''^T=,JiPi^-k) = 0, 1 = l,...,m- 1 

fe=0 

Denoting e~**^i=o •^'^'^^'^^ by 77^, the above constraints are rewritten as follows 

m 

^Pi(xfc)77fc7fc = 0, i=l,...,m-l, 

fe=0 

ni 

k=0 

rn 

T.P^i^k)riklk^ f. (3-23) 

fe=0 
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From invertibility of the matrix Pj/j = Pi{xk) (see Ref. [2]) one can rewrite the Eq. (3-23) as 



rjolo 
Vili 







p-i 



(3-24) 



Vd-ild-i 

V Vdid J \ f ] 

The above equation imphes that rjk^k for k — 0, 1, ...,m are the same as the entries in the 
first column of the matrix P^^ = WP^ multiphed with / and the entries in the last column 
multiplied with /', i.e., the following equations must be satisfied 



Vk^k = 7fee-^*Sr:o^^^K-.) = {WP'),J + {WP').f , ior k = 0, 1, ...,m, 



(3-25) 



with W := diag {^0,^1, ... ,jra)- By using the fact that 7^ and (W^P*)fcm ^^^^ ^ ~ 
0,1,..., m, and so we have -fk = \{WP^)kJ and 7^ = (W^P*)feo- The Eq. (3.26) can be 
rewritten as 

r)k = e-'^^toWi^^) = / + a{k)f (3-26) 

where a{k) is defined as 

{—1 for odd k 
(3-27) 
1 otherwise 

Assuming / = |/|e'^ and /' = |/'|e*^', it should be considered 9' ^ 9 ± ^ then 

e-*Er=o^'^Kx.) ^ e^^d/l ±i<j(A;)|/'|) = e^(^±axctan(^)+2c..). ^^^^ (3_28) 

One should notice that, the Eq. (3.29) can be rewritten as 

1 a(k)\f'\ a(k)\f'\ 

(Jo, Ji, . . . , Jm) = -t[^ + 2co7r ± arctan ( — — — ), 9 + 2ci7r ± arctan ( — — — ), 
* 1/1 I/I 

,^ + 2c^7r±arctan(^^^^)](W^P*) (3-29) 



or 



I/I 



•^fc = -7 + 2c,-^ ± arctan (^i^)](W-P*)^., 
j=o I/I 



(3-30) 
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where Cj for j = 0, 1, . . . , m are integers. The result (3-30) gives an exphcit formula for suitable 
couphng constants so that GHZ state in the final state can be achieved. The final state is as 
the form 

iV'(i)) = /|11...100...0) + /'|00...011...1) (3-31) 

One attempt to provide a computationally feasible and scalable quantification of entanglement 
in multipartite systems was made in Refs. [40,41,42]. For a pure n-qubit state the so-called 
global entanglement is defined as 

-1 N-l 

g(lV')) = 2(1-- 5: Tr[p,^]) (3-32) 

-'^ i=0 

where pi represents the density matrix of ith qubit after tracing out all other qubits. As 
seen from this definition, the global entanglement can be interpreted as the average over the 
(bipartite) entanglements of each qubit with the rest of the system. The global entanglement 
for state in Eq. (3.32) will be 

Qm) = mvf (3-33) 

Also we introduce a simple multiqubit entanglement quantifier based on the idea of bipar- 
tition and the measure negativity (which is two times the absolute value of the sum of the 
negative eigenvalues of the corresponding partially transposed matrix of a state p) [43]. For 
an arbitrary A^-qubit state Psis2-SN > ^ multiqubit entanglement measure can be formulated 

as [44] 

JV 

N 2 

^ = y Qk\N-kiPsis2...SN) (3-34) 

where N is assumed even, otherwise y should be replaced by and gk\N-k{PsiS2-SN) 
entanglement in terms of negativity between two blocks of a bipartition k\N — k of the state 
PsiS2-SN- We can define the following partition-dependent residual entanglements (PREs) 

TT — ^ 

'■'■qi—gm.qm+l---qk\<lk+l---<lngn+l---qN ~ Sqi—qmqm+l—qk\qk+l—qnqn+l—qN 
~^qi—qm\qk+l—qn ~ Qqi...qm\qn+l—qN ~ Sqm+i...qk\qk+i—qn ~ ^qm+l—qk\qn+l—qN (3-35) 
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and 

k N 

qi-qk\gk+i-gN ~ ^qi-qk\qk+i-qN -EE 4.. (3-36) 

i-1 j=k+l 

For the state in Eq.(3.32), we have 

^qi-qmqm+l-qk\qk+l-qnqn+l-qN — ^qi...qk\qk+i-qN ~ 6qi...qk\qk+l-qN ~ ^l/l 1/ I (3-37) 

Another useful entanglement measure was introduced in Refs. [45,46] for n-qubit state 1-0) = 
X^^^o ^ ^iV) with even n, as 

T(^) = 2|x*(a,n)| (3-38) 

where 

X*{a,n)= J2 sgn*{n,i){a2ia(2n-i-iy2i - a2i+ia(^2n-2-2)-2i), (3-39) 

i=0 

( (-1)^« < i < 2"-3 - 1 
sgn*{n,i) = <^ (3-40) 

[ (^-l^Nii)+n 2"-3 < i < 2"-2 - 1 

where, N{i) is the number of the occurrences of 1 in the n-bit binary representation of i as 
in-i---iiio ( in binary representation, i is written as i = i„_i2"~^ -|- ... -|- ii2^ + io2°). For the 
state Eq.(3.32), one can see that 

r(^) = 2\x*{a,n)\ = 2|a2™-ia22™_2-| = 2|//'|. (3-41) 

In order to achieve maximal entanglement [GHZ state), we should have 

I/I = 1/1 = ^ (3-42) 

Then Q(\lp)) = ^qi...qmqm+i---qk\qk+i-qnqn+i-qN ~ ^qi---qk\qk+i-qN ~ '^i'^) ~ ^^ 

In the following we consider the four qubit state (the case m — 2) \tl;{t — 0)) = jllOO) and 
the six qubit state (the case m = 3) \y'j{t = 0)) = 1 111000) in details: From Eq. (2-12), for 
m = 2, the QD parameters are given by 

ai = 2, Q!2 = 0; 001 = 002 = 4, 
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Then by using the recursion relations (B-48) and (B-51), we obtain 



Q^^\x) =x^ -2x-4:, Qaix) = x{x - 4){x + 2), 



so that the stieltjes function is given by 

Q^\x) 



x'^ -2x- A 



Qz{x) x{x-A){x + 2)' 
Then the corresponding spectral distribution is given by 

2 1 

= E7/5(a^ - xi) = -{35(a;) + 5{x - 4) + 28{x + 2)}, 

1=0 " 



which indicates that 



W 



70 ^ 
7i 
\Q 72 y 



1 ^ 



3 
2 



In order to obtain the suitable coupling constants, we need also the eigenvalue matrix P with 
entries Pij = Pi{xj) = --j==Qi[xj). By using the recursion relations (B-48), one can obtain 
Po(^) = 1, Pi{x) = f and P2{x) = -2x- 4), so that 



^ 1 1 ^ 



2 0-1 

1 -1 1 



Then, Eq. (3-30) leads to 



TT 



-i( Jo + 2 Ji + J2) = ^ ± - ± 2co7r, 



TT 



-t{Jo -J2)^e^-± 2ci7r, 



TT 



-t( Jo - Ji + J2) = ^ ± - ± 2c27r. 



Now, by considering co = ci = C2 = we obtain 



_ e _ _ TT 

Jo — —J: Jl — 0, J2 — 
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Also by considering cq = 0, ci = C2 = 1 the coupling constants will be 



, J^ — it , Jo — it 

3i ' 3i' m 



and by considering cq — l,ci — C2 — 

3^' -^^-^S^' -^'-^Ut 
From Eq. (2-12), for m = 3, the QD parameters are given by 

CKi =4, q;2 = 4, as = 0; c^i = o's = 9, uj2 — 16 

Then by using the recursion relations (B-48) and (B-51), we obtain 

= -8x^ -9x + 36, Q4x) = (x^ - 9)(x - 9)(x + 1), 



so that the stieltjes function is given by 



x^ - 8a;^ - 9a; + 36 



Qi{x) {x"^ - 9){x - 9){x + I) 
Then the corresponding spectral distribution is given by 

^ 1 

H[x) = Y,-fid{x - xi) = —{S{x - 9) + 56{x - 3) + 9S{x + 1) + 56{x + 3)}, 
1=0 20 



which indicates that 



W 



70 ^ 

7i 

72 

\^ 73 y 



1 

20 



'^ 1 ^ 



5 
9 
5 



By using the recursion relations (B-48), one can obtain Po{x) — 1, Pi(x) = |, P2{x) — 
^(x^ - 4a; - 9) and Psix) = ^{x^ - 8x^ - 9x + 36), so that 



1 1 

3 1 



3-1-1 1 
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Then, Eq. (3-30) gives 



t{Jo + 3Ji + 3J2 + J3) 



--0±-± 2co7r, 

Ot^T 2ci7r, 
= 0±j± 2c27r, 
Ot^± 2c37r. 



t{Jo + J1-J2- J3) 



t{Jo 




t{Jo -J1 + J2- J3) 



Again, by considering cq 



Ci — C2 — C3 — we obtain 



Jo — ~T) Ji — J2 — 0, J3 — — -. 



4 Conclusion 



A 2m-qubit initial state was prepared to evolve under a particular spin Hamiltonian, which 
could be written in terms of the adjacency matrix of the Johnson graph J(2m,m). By using 
spectral analysis methods and employing algebraic structures of the Johnson networks, such as 
distance-regularity and stratification, a method for finding a suitable set of coupling constants 
in the Hamiltonians associated with the networks was given so that in the final state, the 
maximal entanglement of the form GHZ state, could be generated. In this work we imposed a 
constraint so that all amplitudes in the final state were equal to zero except to two amplitudes 
corresponding to the first and the final strata (any pair of antinodes of the network), where 
for J(2m, m) these strata contain only one vertex, then GHZ state was generated. We hope 
to generalize this method to arbitrary Johnson networks J(n, m) and other various graphs, in 
order to investigate the entanglement of such systems by using some multipartite entanglement 
measures. 
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Appendix 

A Stratification technique 

In this section, we recall the notion of stratification for a given graph F. To this end, let d{x, y) 
be the length of the shortest walk connecting x and y for x ^ y. By definition d{x, x) = 
for all X E V. The graph becomes a metric space with the distance function d. Note that 
d{x, y) = 1 ii and only ii x ^ y. We fix a vertex o G as an origin of the graph, called the 
reference vertex. Then, the graph F is stratified into a disjoint union of strata (with respect 
to the reference vertex o) as 

oo 

V^\Jri{o), ri{o):^{aeV :d{a,o)^i} (A-43) 

1=0 

Note that Fi(o) = may occur for some i > 1. In that case we have Fj(o) = Fi+i(o) = ... = 0. 
With each stratum Fj(o) we associate a unit vector in 1^{V) defined by 

10.) = ^ E (A-44) 

where, Ki — |Fi(o)| is called the i-th valency of the graph {Ki := \{j : d{o,j) — i}\ = |Fi(o)|). 

One should notice that, for distance regular graphs, the above stratification is independent 
of the choice of reference vertex and the vectors = 0,1, ...,d — 1 form an orthonormal 

basis for the so called Krylov subspace Kd{\(f)o),A) defined as 

Xd(|0o),A) =span{|(/)o),A|0o),---,A^-^|0o)}. (A-45) 
Then it can be shown that [25], the orthonormal basis are written as 

|0,) = P,(>l)|0o), (A-46) 

where Pj = ao + aiA+ ... + aiA^ is a polynomial of degree i in indeterminate A (for more details 
see for example [25,26]). 
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B Spectral distribution associated with the graphs 

In this section we recall some facts about spectral techniques used in the paper. From or- 
thonormality of the unit vectors given in Eq.(A-44) (with |0o) as unit vector assigned to 
the reference node) we have 

Sij = {<i>i\4>j) = I P^{x)P^{x)^l{dx). (B-47) 

By rescaling as = y/uJ^TTTuJkPk, the spectral distribution fi under question will be 
characterized by the property of orthonormal polynomials {Qk} defined recurrently by 

Qo{x)^l, Qi{x)^x, 

xQk{x) = Qk+i{x) + OikQk{x) + UkQk-i{x), k>l. (B-48) 
The parameters ak and cuk appearing in (B-48) are defined by 

ao = 0, ak = K-hk- Ck, ujk = = h-iCk, k = l,...,d, (B-49) 

where, k = ki is the degree of the networks and biS and q's are the corresponding intersection 
numbers. Following Refs. [34], we will refer to the parameters ak and ouk as QD (Quan- 
tum Decomposition) parameters (see Refs. [26,27,28,34] for more details). If such a spectral 
distribution is unique, the spectral distribution /i is determined by the identity 



'^^^^^ Irx-v x-ap- ^ ^ Qd+i{x) (B 50) 

x — an — " 

X — a3 

where, xi are the roots of the polynomial Qd+i{x). G^{x) is called the Stieltjes/Hilbert trans- 
form of spectral distribution /j, and polynomials {Q^ ^} are defined recurrently as 

Q^^\x) = l, Qf\x)=x-a,, 

xQ^,}\x) = Q^l,{x) + ak+iQ^\x)+ujk+iQ^l^{x), k > 1, (B-51) 
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respectively. The coefficients 7/ appearing in (B-50) are calculated as 

7/ := [(x - xi)Gf,{x)] (B-52) 

Now let G^{z) is known, then the spectral distribution /j, can be determined in terms of 
xi, 1 — 1,2, ... and Gauss quadrature constants 7;, 1 — 1,2, ... as 

d 

t^^^liK^-^i) (B-53) 

l=Q 

(for more details see Refs. [35,36,37,38]). 
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